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Abstract 

'We report a detailed magnetic study of a new type of self-organized nanowires disclosed briefly pre- 
viously [B. Borca et ai, Appl. Phys. Lett. 90, 142507 (2007)]. The templates, prepared on sapphire 
wafers in a kinetically-limited regime, consist of uniaxially-grooved W(llO) surfaces, with a lateral 
period here tuned to 15 nm. Fe deposition leads to the formation of (110) 7nm-wide wires located at 
|the bottom of the grooves. The effect of capping layers (Mo, Pd, Au, Al) and underlayers (Mo, W) 
X)n the magnetic anisotropy of the wires was studied. Significant discrepancies with figures known for 
thin fiat films are evidenced and discussed in terms of step anisotropy and strain-dependent surface 
anisotropy. Demagnetizing coefficients of cylinders with a triangular isosceles cross-section have also 
been calculated, to estimate the contribution of dipolar anisotropy. Finally, the dependence of mag- 
petic anisotropy with the interface element was used to tune the blocking temperature of the wires, 
here from 50 K to 200 K. 
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1. Introduction 

Bottom-up processes are subject to an increas- 
ing interest for the fabrication of nanostructures 
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at surfaces. These processes are inexpensive alter- 
natives to the ever-rising cost of lithography with 
the decreasing pitch. Self-organized dots or wires 
i.e. arranged as a regular array on surfaces, are of 
particular interest because a narrow size distribu- 
tion may arise as a consequence of the order [1]. 
Such studies are motivated by the demand for 
an ever higher density of hard-disk drives, for 
which goal reducing the dispersion of grain size 
and magnetic properties of the media is a require- 
ment. Beyond this technological motivation self- 
organized systems contribute to the fundamen- 
tal understanding of magnetism in low dimen- 
sions, based on chemical[2] or epitaxial routes[3]. 
Nanostructures from the micron size [4, 5, 6] down 
to the atomic size [1, 7, 8, 9] have been studied, 
for studying in low-dimensions e.g. the micromag- 
netics of magnetic domains and domain walls, and 
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magnetic moments and anisotropy, respectively. 

We focus here on the fabrication and mag- 
netic properties of self-organized arrays of epitax- 
ial nanowires. This system addresses the tune- 
ability in terms of period and magnetic properties, 
which are one of the prerequisites for such systems 
to be relevant for applications. Magnetic wires at 
surfaces are often achieved by step-decoration of 
vicinal surfaces[10, 11, 12]. In this approach a 
crystal has to be polished with a specific miscut 
to achieve the desired period. Templates result- 
ing from kinetic effects are potentially more versa- 
tile as the period can be changed with parameters 
like temperature. Ion etching under grazing inci- 
dence may be used to create ripples and wires on 
surfaces[13]. However a significant control of the 
period has not been demonstrated for magnetic 
materials so far. Here we report on the advance- 
ment of a recent alternative approach [14, 15], 
based on the growth of a body-centered-cubic 
W(llO) template on nominally-flat sapphire, in- 
stead of etching. Fe wires are then grown at the 
bottom of the trenches in a layer-by-layer fash- 
ion. Using the dependance of surface magnetic 
anisotropy with the interfacial element, we could 
tune the magnitude of the magnetic anisotropy of 
the wires grown on such templates using suitably- 
chosen capping layers and underlayers. As a con- 
sequence this gives us a mean to tune the mean 
blocking temperature of the wires, here from 50 K 
to 200 K. 

2. Experimental procedures 

The samples were grown in a set of ultra- 
high vacuum chambers, using pulsed-laser depo- 
sition with a Nd-YAG laser operated at 532 nm. 
The metallic films and nanostructures are grown 
starting from commercial sapphire single crystals. 
The typical fiuence is 1 — 3 J.cm~^.pulse~^ de- 
pending on the element, yielding a rate of de- 
position around lA/min. The deposition cham- 
ber is equipped with a quartz microbalance, 
sample heating and a translating mask for the 
fabrication of wedge-shaped samples. We use 
a Riber lOkeV Reflection High Energy Elec- 
tron Diffraction (RHEED) setup with a CCD 
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Figure 1: (a,b) 100 nm x 100 nm STM images of (a) the 
W(llO) template (b) w Inm-thick wires of Fe(0 = 
2.5 AL) /W. (c) Cross-sections along the lines shown in (a)- 
gray and in (b)-black 



camera synchronized with laser shots, which 
thus permits operation during deposition. An 
Omicron room-temperature Scanning Tunneling 
Microscope (STM), used in the constant cur- 
rent mode, and an Auger Electron Spectrome- 
ter (AES) with a MAC-II analyzer are located in 
connected chambers. A quantitative analysis was 
performed based on the energy range 100 — 800 eV 
and an incident beam energy of 5keV, and fo- 
cusing on the NOO and LMM transitions of W 
and Fe, respectively centered around 180 eV and 
650 eV. Each measured spectrum was first de- 
rived, then smoothed and normalized by the back- 
ground, representing the signal of secondary elec- 
trons detected in the energy range where no char- 
acteristic peaks were expected. A detailed de- 
scription of the chambers and growth procedures 
can be found in [16]. The magnetic measure- 
ments were performed ex situ with a Quantum 
Design MPMS Superconducting QUantum Inter- 
ference Device (SQUID) magnetometer. 

3. Nanowires preparation 

The ordered arrays of nanowires are obtained 
in several steps, described below. 

The first step is the preparation of a smooth 
metallic buffer layer on sapphire (1120) wafers. A 




Figure 2: a) Schematic view of the sample prepared for the quantitative AES analysis, starting from two areas with a 
flat and a grooved W(llO) surface. The two ends provide thick layers of W and Fe, connected by a continuous wedge 
of Fe b) Reference spectra of Fe (i?Fe) and of W (Rw) respectively c) Spectra at position x along the wedge and the 
fit as a linear combination of the reference spectra d) Variation of the fitted intensities cw and cpc of the W NOO and 
Fe LMM peaks over the wedge, for both the flat and the grooved surfaces e) Schematic representation of the expected 
curves for a Frank van der Merwe, Stranski-Krastanov and Volmer Weber growth. 



seed layer of Mo (nominal thickness 6 = 1 nm) 
followed by W (9 ~ 10 nm) are deposited at 
room temperature (RT) followed by annealing 
at 800 °C. This yields a smooth (110) surface, 
with atomically-fiat terraces of mean width k, 
200nm[16, 17]. 

The second step consists of the preparation of 
the non-magnetic self-organized template. This 
step was inspired by reports of a kinetic uni- 
axial roughening of films of the bcc elements 
Fe(110)[18] and W(110)[19] during homoepitaxy 
at moderate temperatures. This roughening was 
explained by anisotropic diffusion along atomic 
steps and the occurrence of an Ehrlich-Schwoebel 
barrier. Based on these clues we optimized the 
deposition parameters to produce self-organized 
arrays of trenches aligned along the [001] direc- 
tion of the W(llO) surface (Fig. la). Here we focus 
on arrays deposited at 550 °C for the first atomic 
layers, and progressively reduced to 150 °C during 
growth. A stable period is then reached within a 
few nanometers of deposition. The period and 
depth of the trenches were determined by STM 
to be ~ 15 nm and ~ 2.5 nm for these deposition 
parameters, respectively. The trenches have well- 
defined facets of {210} type, revealed by STM and 
RHEED (tilted streaks) with a constant angle of 
9 ^18° with respect to the mean surface[20, 14]. 
The occurrence of facets with a well-defined an- 
gle minimizes the effect of the fluctuations of the 
width of the trenches on the shape of the trenches, 
which thus displays essentially no distribution. 



This is desirable as the distribution of magnetic 
anisotropy energy may increase tremendously for 
nanosized systems in relation with the distribu- 
tion of structural environments [21]. 

As an optional third step one atomic layer of 
Mo can be deposited at 150 °C on this template, 
allowing a control of the chemical nature of its 
free surface without affecting the morphology of 
the template nor the angle of the facets. 

As the next step these self-organized templates 
are used for fabrication of magnetic wires. To 
this end we deposit Fe at 150°C, a temperature 
under which Fe grows essentially layer-by-layer on 
nominaUy-Bat W(110)[18]. STM reveals a level- 
ing of the bottom of the trenches, suggesting their 
progressive filling and thus presumably yielding 
Fe wires with a triangular cross-section (Fig.lb- 
c). This conclusion is supported by the quanti- 
tative analysis of AES spectra. To avoid errors 
in thickness resulting from fluctuations in depo- 
sition rates from one sample to another, all mea- 
surements were performed on one single wedged 
sample as depicted in Fig. 2a. Two ends of the 
sample consist of a 3nm thick Fe layer and of 
a 9 nm thick W layer respectively, which provide 
reference spectra Rpc and -Rw for pure elements as 
well as a measure of the beam intensity (Fig. 2b). 
In-between these two areas Fe is deposited with 
a uniformly-varying thickness, both on a grooved 
and on a fiat W(llO) surface. Spectra at position 
X along the wedge were fitted linear combi- 
nation of the reference spectra (Fig. 2c) : 



Figure 3: Magnetic measurements of Fe/W(110) wires with a capping layer of Mo; upper-part without under-layer, 
lower-part with Mo under-layer, a) Magnetization curves at 10 K in-plane along [001] and [110]. b) Magnetization 
curves in-planc along [001] at different temperature, c) Dependance of coercivity with temperature, along with a \/T 
fit. d) Zcro-ficld-coolcd (ZFC)/ficld-coolcd (FC) curves, measured with a low magnetic field (0.01 T) applied along the 
in-planc [001] direction. 



l{x) = Cw(a;)-Rw + CFe(x)i?Fe (1) 

where cy/{x) and cfc{x) are the intensity of W 
and Fe. Fig.2d shows the variation of cyj{x) and 
cpeix), on the flat and the grooved parts, with the 
Fe thickness x varied in the range 0-1 AL (Atomic 
layers). First notice that the mean free path of 
electrons at the MNN Fe peaks is higher than 
that at the NOO W peaks, explaining a faster 
decrease of cw than the increase signal of cpc- 
Then, is it clear that the variation of the cw(a^) 
and CFe{x) is faster for the flat surface, where 
Fe is known to grow layer-by-layer (Frank van 
der Merwe mode), than for the grooved template. 
This behavior, reminiscent of the Volmer- Weber 
growth mode (Fig.2e) with three-dimensional is- 
lands on a bare substrate [22], may suggest that 
Fe does not or not fully wet the microfacets of 
the W template and flows from the beginning to 
the bottom of the grooves to form wires. 

For the magnetic measurements reported here 
the growth of the wires has been stopped at 
= 2.5 AL, i.e. before the percolation thresh- 
old. The mean full width and mean full height 
are then expected to be 7nm and 1.2 nm, respec- 
tively. Notice that owing to the expected trian- 
gular shape of the wires, the height averaged over 
a wire equals half its full height. 

The final step is the deposition at room tem- 
perature of various capping layers (Au, Al, Mo, 



Pd), with the aim of modifying the magnetic 
anisotropy. The thickness of this layer is typically 
=4-5 nm, to act also protection against 
oxydation or other adsorbate-induced contamina- 
tion, except for Pd where only an interface layer 
of thickness =4 AL is deposited to set the mag- 
netic anisotropy, capped by 4 nm of Au as a pro- 
tection against contamination. 

4. Magnetic measurements and analysis 

The mean period, full width and full height of 
the wires investigated are 15 nm, 7nm and 1.2 nm, 
respectively. SQUID hysteresis loops were per- 
formed on assemblies of such wires along three 
directions: along the direction normal to the 
plane, and along two in-plane directions: paral- 
lel to Fe[001] (along the wires) and parallel to 
Fe[110] (across the wires). 

We studied the influence on hysteresis loops of 
the chemical nature of both the top and of the 
bottom interfaces. For the top interface the ele- 
ments Pd, Au, Al and Mo have been used (Fig.4). 
For the bottom interface only W and Mo have 
been used so as to keep the growth process of 
the wires unchanged (above room temperature 
Fe is liable to surface alloying or intermixing 
when deposited on Al, Pd or Au). In all cases 
the easy axis of magnetization essentially remains 
along [001], i.e. along the length of the wires. 
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Figure 4: Magnetization curves at 10 K in-plane along 
[001] and [110], and perpendicular to the plane along [110] 
for Fe wires of type: a) Pd/Fe/W; b) Pd/Fe/Mo/W; c) 
Au/Fe/W; d) Al/Fe/W. 



The nature of both interfaces affects the coer- 
civity along this direction, and the anisotropy 
fields along the two harder directions, in-plane 
[1 — 10] and out-of-plane [110]. The experimental 
density of anisotropy energy along a hard direc- 
tion was calculated from the hysteresis loop as 
-E'exp = A^o /o ' HdM. To minimize systematic 
errors related to any small opening of the loop 
(caused e.g., by the misalignment of the sample 
in the Squid) we processed the anhysteretic curve 
(sum of the two branches Mup and Mdown) instead 
of one single branch. The experimental values 
thus determined of both in-plane -ft'tot.p and out- 
of-plane Ktot,e densities of anisotropy energy are 
summarized in Tab.l. In the case of Pd/Fe/W 
hysteresis is observed in both in-plane directions, 
which reveals that the anisotropy is weak and we 
are on the verge of a spin reorientation transi- 
tion in-the-plane. In this case a significant error 
probably arises from the crude procedure of using 
the calculated anhysteretic curve for fitting the 
anisotropy. 

Let us analyze these values, with a final discus- 
sion on the magnitude of surface anisotropy with 
respect to values known for fiat thin films. The 
total anisotropy originates as the sum of several 
contributions. Some of them can be estimated, 
while others are entangled and their distinction is 
mainly a view of the mind. We use the following 
indices to label anisotropy constants and demag- 
netizing factors: no index for the total energy, d 
for dipolar energy (both expressed in J/m^), s for 
surface (expressed in J/m^), p for in-plane (all 
with the convention that positive figures mean 
easy axis along [001]), 9 for out-of-plane (with the 
convention that negative means easy axis perpen- 
dicular to the plane). Notice that unlike some of 
the existing works [26, 25] we use the same sign 
convention for all anisotropics (total, volume and 
surface) . 

The first type of anisotropy is the bulk fourth- 
order anisotropy of Fe, with Ki = 4.58 x 10^ J/m'^. 
When projected along a plane of the (< 110 > 
, < 001 >) family this term yields a second or- 
der contribution Ki and a fourth order contribu- 
tion — {3/A)Ki [25], whose sum is negligible before 
most figures found experimentally here. 



Table 1 : Experimental in-plane (labeled p) magnetic anisotropy energy density for various capping or imdcrlayers, derived 
at 10 K from the hysteresis loops of Fig. 3 and Fig. 4. The numbers shown in brackets are theoretical ones, computed with 
values known for flat (110) thin films [23, 24, 25]. Concerning the differences only those numbers are displayed, which 
allow a comparison with a change of one interface only with Mo/Fe/W or Mo/Fc/Mo/W. Out-of-plane anisotropy has 
not been studied extensively, so that only the available flgures are indicated. 

Interfaces i^tot.p (10^ J/m^) AKs,p (mJ/m^) with AJsTs.p (mJ/m^) with A'tot.e (10^ J/m^) 

Mo/Fe/W Mo/Fe/Mo/W 

Mo/Fe/W 3.6 ±0.6 — -0.03 ± 0.04(-0.56) 

Mo/Fe/Mo/W 4.2 ± 0.4 0.03 ± 0.04(0.56) — 

Pd/Fe/Mo/W 1.2 ±0.2 -0.12 ± 0.03(0.18) -0.15 ± 0.02(-0.38) 16.3 ± 0.06 

Pd/Fe/W 0.35 ±0.15 -0.16 ± 0.03(-0.38) — 7.1 ±0.1 

Au/Fe/W 1.7 ±0.3 -0.10 ± 0.03(0.21) — 7.4 ±0.1 

Al/Fe/W 0.9 ±0.05 -0.14 ± 0.03(-0.24) — 



The second type is dipolar energy, arising from 
both the self energy of a wire, and from the in- 
teractions within the array. With the idealized 
picture of equally-spaced infinite cylinders with 
a triangular isosceles base, let us call A'w.p and 
A'w.e the transverse and perpendicular demagne- 
tization coefficients of an isolated wire, and j?Va,p 
and A^A,0 the coefficients for the entire array (also 
including the self-dipolar energy of each wire). 
A^v,p ± N-w^g = 1 and A^a,p ± Na,9 = 1 because 
the wires are assumed to be infinitely long in the 
[001] direction. Based on analytical calculations 
one finds iVw,p — 0.147 (see Annex 1) with the 
present geometry. The effect of the array is only 
minor and one finds, again analytically(Annex 2): 
Na,p — 0.142. Based on bulk Fe magnetization 
Ms = 1.73 X lO^A/m demagnetizing energy den- 
sities are: -ft'd.w.p = 2.76 x 10^ J/m^ and -ft'd.w.e = 
1.60 X lO^J/m^, and i^d,A,P = 2.67 x lO^'j/m^ 
and Kd,A,e = 1-61 x lO^J/m^. Thus dipolar en- 
ergy arises primarily from self-dipolar energy of 
individual wires, and is in magnitude and sign a 
major contribution to the value measured exper- 
imentally. 

Another contribution to magnetic anisotropy is 
of magneto-elastic origin, both linear and non- 
linear [27]. However unlike thin films investigated 
so far, stress is both in-plane and out-of plane 
owing to the occurrence of many steps, imply- 
ing a significant shear because W and Fe have 
a 10 % lattice mismatch. Therefore no figure or 
even sign can be reliably predicted, as no reference 



figure are available for the non-linear corrections 
to magneto-elasticity. 

All three types of anisotropics mentioned above 
are of volume type and to first approximation 
should be identical for all systems reported here, 
of identical thickness. With this respect, notice 
that the 1 AL-thick Mo underlayer is pseudomor- 
phic with W so that no difference in magneto- 
elastic contribution is expected with or with- 
out the underlayer. Then there remains surface 
anisotropics, which should thus be responsible for 
any change found between the different samples. 
Whereas absolute values of interfacial anisotropy 
can not be derived for wires, differences of inter- 
facial energy between wires where only the top 
interface is changed, can be derived from these 
experiments (Table 1). These differences are dis- 
cussed in the following. 

Let us first notice the consistency within the 
present various experimental data. Experimen- 
tal figures are available for the change from Mo 
to Pd capping interface for both W and Mo 
underlayer. The resulting changes of surface 
anisotropy are in very good agreement one with 
another, — 0.16mJ/m^ and — 0.15mJ/m^, respec- 
tively. This noticed, let us discuss the results with 
respect to the literature. The change of magnetic 
anisotropy upon changing the interface, computed 
from published data concerning thin films, are 
shown in brackets in Table 1. The experimen- 
tal (wires) and computed (films) figures concern- 
ing the change of the bottom interface (first three 
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Figure 5: ZFC/FC curves measured with a low magnetic 
field (0.01 T) applied along the [001] direction for: a) 
Pd/Fc/Mo/W; b) Al/Fe/W. 



lines of Table 1) are in clear disagreement. This 
may point at the role of the atomic steps, which 
are known to induce an extra anisotropy energy 
per atom with respect to a flat surface[28, 29]. 
The step-edge anisotropy for Fe/W(110) with 
steps oriented along the [001] direction was deter- 
mined to be i^step,p ~ -1 X 10"^3 J/m[28]. With 
one step every other 6.7 A for {210} facets, this 
yields 0.2 mJ/m^ as an order of magnitude for the 
extra contribution to the interfacial anisotropy 
energy. This figure is of the same order of magni- 
tude as the discrepancy between our experimen- 
tal data for the stepped underlayer interface and 
the data computed from fiat thin films. How- 
ever no experimental data is available for vici- 
nal Mo surfaces so that the discussion cannot be 
pushed forward. The figures for the change of 
anisotropy energy related to cappings layers are 
also in disagreement (see especially the case of 
Au) , although no change of step anisotropy is ex- 
pected in this case. This may hint to a depen- 
dance of surface anisotropy with strain, here both 
in-plane and out-of-plane, once claimed[30] and 
since revisited in terms of non-linear magnetoe- 
lastic coupling coefficients in the simple situation 
of uniform strain[31]. 

Concerning the out-of-plane anisotropy, experi- 
mental figures are available for Pd and Au capping 
layers. A very small difference in the experimen- 
tal values of the magnetic anisotropy is obtained 
for Pd/Fe/W and Au/Fe/W systems (Table 1). 



This feature is in good agreement with that of 
the thin films, contrary to the case of in-plane 
anisotropy (discussed previously). With the same 
reassessment as for the in-plane anisotropy and 
considering the known data from thin films for 
K^fi{kn/¥e)= -0.72 mJ/m^ and K,^g{Vd/¥e) = 
— 0.66mJ/m^ [25] it is figured out that magneto- 
elastic terms and deformations favor the out-of- 
plane configuration with an anisotropy energy of 
2.2 X 10^ J/m'^. The observed behavior points at 
the dependence of the surface anisotropy with the 
strain. In this sense, the STM and RHEED mea- 
surements (not presented here), because no dis- 
locations where observed, suggest that the wires 
are pseudomorphic up to thicknesses around 1 nm 
higher than that of the thin films (< 2 ALs). This 
has probably a major impact on magneto-elastic 
anisotropy. 

We finally discuss some aspects of coercivity, 
which unlike anisotropy is an extrinsic phenomena 
and as such may yield information on e.g., the mi- 
crostructure of a system. We focus on the case of 
Mo-capped Fe/W wires. At low temperature the 
hysteresis loops measured with the external field 
H/ /[OOl] have a rather square shape, and at rema- 
nence full saturation is still observed (Figure 3b). 
The coercivity continuously decreases with tem- 
perature while Mg remains essentially unchanged 
suggesting a superparamagnetic behavior when 
coercivity vanishes. The temperature-dependent 
coercivity HdT) is plotted in Figure 3c. The de- 
crease of Hc{T) reasonably follows a \/T law, as 
expected for a uniaxial anisotropy K within the 
framework of the coherent rotation model, assum- 
ing in a first approximation that K does not de- 
pend on temperature. Under these two hypothe- 
ses the energy barrier preventing the reversal of 
magnetization varies like 



AEoc [1- H/H^{T = OK)f 



(2) 



and assuming a thermally activated magnetiza- 
tion switching with a switching rate following an 
Arrhenius law with a time constant 



T = Toexp{AE/kBT), 



(3) 



one derives: 



Table 2: The blocking temperature Tb determined from equation (3) and (4) for Hc{Tb) = and in brackets from ZFC 
curves respectively, the nucleation volume and the corresponding nuclcation length for various capping or underlayers, 
and the domain wall width 



Interfaces 


Tb(K) 


Vn(nm^) 


/(nm) 


wall width (nm) 


Mo/Fe/W 


110 (100) 


75 (68) 


18 (16) 


23 


Mo/Fe/Mo 


200 (190) 


165 (157) 


39 (37) 


22 


Pd/Fe/Mo 


(205) 


(590) 


(140) 


41 


Al/Fe/W 


(50) 


(190) 


(45) 


47 



2K 



1 - 



ln(r/ro)A:BT 
KV 



(4) 



The blocking temperature Tb determined as 
Hc{Tb) = is ~ 110 K. This agrees well with 
the half-way-up of the zero-field cooling remag- 
netization curve (found at 100 K), which should 
represent the mean value of the blocking tem- 
perature. Using the relationship KV^^ — 2bk-QT-Q 
where Vn is the nucleation volume, inferred from 
an attempt frequency tq ~ 10"^*^ s and a mea- 
surement duration r ^ 10 s, we conclude that 
~ 75 nm^, with the experimental anisotropy 
energy and the Tb determined from the Hc{T) 
variation. Based on the estimated cross-section 
area of the wires (4.2 nm^), is converted into 
a length of ~18nm. This length is comparable 
to the domain wall width tc^J A/ K=2'imiv. This 
suggests that the Mo/Fe/W wires behave essen- 
tially like individual wires. A further argument 
enforcing this conclusion is the progressive remag- 
netization upon heating during the FC measure- 
ments. Both a much larger and a sharp tran- 
sition in the rising ZFC curve occur in the case of 
strongly coupled wires [32]. 

In the view of this analysis of coercivity the de- 
pendence of the magnetic anisotropy on the un- 
derlayer and capping layer can be used to tune the 
blocking temperature. This is illustrated on Fig- 
ure 5 were systems with high and low anisotropy 
yield a change of mean blocking temperature by 
a factor of four, with an identical volume. 

Table 2 summarizes the values of Tb deter- 
mined from Hc{Tb) = variation and/or ZFC- 
FC curves, and the resulting nucleation volume 
for different available systems. The values of 



the length obtained from the conversion of the V^, 
remain comparable to the domain wall width for 
almost all the cases. Thus the wires are essentially 
isolated and have an individual magnetic behav- 
ior. An exception occurs for Pd/Fe/Mo/ W, when 
the corresponding nucleation length ^ 140 nm 
is much higher than the domain wall width. 
The coupling between the wires probably arises 
through the Pd layer, due to its high induced 
polarization, as observed in magnetic multilayers 
[33, 34]. 



5. Conclusion 

Arrays of epitaxial Fe(llO) nanowires were ob- 
tained by deposition on kinetically self-organized 
trenched surfaces of W, with a period of 15 nm. 
Owing to a self-limiting effect the nano-facets of 
the trenches are of type {210} with very little an- 
gular spread, which is suitable for reaching a low 
distribution of physical properties such as mag- 
netic anisotropy. Several underlayer (Mo, W) and 
capping layer (Mo, Ag, Pd, Al) elements were 
used and the differences of in-plane magnetic sur- 
face energy was measured. A significant disagree- 
ment was found with figures computed from val- 
ues determined on flat thin films, based on litera- 
ture data. This points at the change of magnetic 
anisotropy owing to steps, already established, 
and possibly the dependance of surface anisotropy 
with strain, which is large in such wires, both 
of tensile and shear type. In the course of the 
discussion of magnetic anisotropy, demagnetizing 
factors for infinite cylinders with an isosceles basis 
were calculated analytically, as well as the impact 
of the ordering as an array on these coefficients. 
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Annex 1: demagnetizing coefficients of a 
cylinder with triangular cross-section 

In this section we calculate the demagnetizing 
coefficients and Ny of an infinite cylinder with 
a triangular cross-section, such as + Ny = 
1^. We do not consider a scalene triangle, how- 
ever restrict the calculation to an isosceles tri- 
angle. Although demagnetizing fields are non- 
homogeneous and the self-consistency of the as- 
sumption of uniform magnetization can rigorously 
not be achieved under any arbitrary external field 
for volumes with surfaces not of second order, de- 
magnetizing coefficients can still be defined via 
the dipolar energy of the hypothetically uniformly 
magnetized state [35]. These coefficients are still 
related to the magnetization energy computed 
from hysteresis loops of infinitely-soft bodies[36]. 
We describe the base isosceles triangle by the 
length i of the two equal sides, and the angle 
a e [0 — 7r/2] between any of the two equal sides 
and the base of the triangle (see inset of Figure 6). 
In the following we use the shortcuts c = cos(a) 
and s = sin (a). 

We evaluate the demagnetizing energy in the 
framework of charge and potential: = 
(/io/2) J cr(f)dS, with H = — grad and a = Mg.n 
(n is the outward normal to surfaces). The poten- 
tial created at the origin by a horizontal linear 
segment of ends i = 1,2 characterized by their or- 
dinate y, abscissas Xj, radii rj and angles a,, and 
bearing a uniform charge a, is 



^With the notation of the experimental part: A^^: = 
A^^vv^p and Ny = N-^^g 




n/8 n/4 3n/8 n/2 

Angle a 



Figure 6: Variation of the demagnetizing coefficients 
N^. (o) and Ny (•) with the base angle a. The geome- 
try of the triangle is depicted in the inset. 



0= ((T/27r) [x2logr2 - Xi logri 

?/(«!- 02) + Xi - X2] (5) 

Assuming uniform magnetization along x only 
the two equal sides of the triangle display mag- 
netic charges, with a linear density a = sMg. 
Making use of Eq. (5), we find after some alge- 
bra: 

= sH^c{c\n2c + sa). (6) 

Upon normalizing by (1/2) fi^M^ and by the 
area i'^cs of the triangle, we finally get the de- 
magnetizing coefficients 

= 2s(cln2c + sa)/7r (7) 
Ny = l-2s(cln2c + sa)/7r. (8) 

Figure 6 shows the variation of N^ and Ny with 
a. We may notice the values expected for limiting 
cases: N^^O) = (thin-film with in-plane magne- 
tization), A^2;(7r/3) = 1/2 (equilateral triangle, i.e. 
with vanishing two-fold anisotropy), Nx{7t/2) = 1 
(thin-film with perpendicular magnetization). 

One easily finds the expansions for fiat triangles 
(a 0): 



41n2t 

vr w 
~ 0.88p 



(10) 
(11) 



and for acute triangles 



a 



vr/2): 



TT 



— a 1 In 



vr 

1 , u; 
m — 

vr t t 

-0.32g Ing 



vr 



— a 



(12) 

(13) 
(14) 



where t = is and w = 2ic are the full thickness 
(along y) and full width (along x) of the trian- 
gle, respectively, and p = t/iy ^ for a — > 
and q = w/t Ofora vr/2 are the trian- 
gle full aspect ratios. These expansions might be 
compared to those of infinite cylinders with an el- 
liptical cross-section with axes lengths w and t: 
Nx ~ p and Ny q for a — s> and a — > 7r/2, 
respectively [37]. This shows that demagnetizing 
factors of infinite cylinders with an isosceles tri- 
angular basis are reasonably well described by an 
ellipsoidal basis with identical encompassing rect- 
angle for a flat triangle, while it is not the case 
for an acute triangle. 

Annex 2: demagnetizing coefRcients of a 
ID array of dipolar lines 

In this section we estimate the corrections that 
shall be made to the demagnetizing coefficients of 
an isolated linear object such as the cylinders with 
triangular cross-section considered above, when 
an infinite one-dimensional array of them is now 
considered. We shall name x the direction of the 
array, with period L (Figure 7). The correction 
will be calculated under the hypothesis that the 
cylinders are sufficiently apart one from another 
so that they can be approximated with linear 
dipoles, i.e. localized magnetic dipoles in a 2D 
space. Here yU. = sM^ where s is the of the 
cross-section of the cylinder, thus s = wt/2 here. 



Figure 7: Geometry and naming for the calculation of the 
dipolar energy in an infinite chain of cylinders with an 
isosceles triangular cross-section. 



Notice that in a 2D space the unit of magnetic 
dipoles is thus A.m. 

It is readily shown that the stray field arising 
from a dipole fj. in a 2D space is 



Hd(r) 



1 



2vrr2 



(15) 



We proceed to calculation based on dipoles di- 
rected along the period of the array, so as to es- 
timate the associated decrease of A''^: \x = /ii. 
Noticing that ^^i(l/''^^) = tt^/S, the stray field 
created at the origin by all other dipolar lines 
reads (/ivr/6L^)i. The correction to is thus: 



vrs 
"6L2 



(16) 



In the case of cylinders with an isosceles cross- 
section, we have AAT^, = — (vr/12)(u;t/L^). 
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